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DRAFT 


ABSTRACT 

We present molecular line imaging observations of three massive molecular outflow 
sources, G333.6-0.2, G333.1-0.4, and G332.8-0.5, all of which also show evidence for 
infall, within the G333 giant molecular cloud (GMC). All three are within a beam 
size (36 arcseconds) of IRAS sources, 1.2-mm dust clumps, various masing species and 
radio continuum-detected Hii regions and hence are associated with high-mass star 
formation. We present the molecular line data and derive the physical properties of 
the outflows including the mass, kinematics, and energetics and discuss the inferred 
characteristics of their driving sources. Outflow masses are of 10 to 40 MQin each 
lobe, with core masses of order 10^ Mq. outflow size scales are a few tenth of a parsec, 
timescales are of several xlO'’’ years, mass loss rates a few xl0“^ MQ/yr. We also hnd 
the cores are turbulent and highly supersonic. 

Key words: stars: formation - ISM: clouds - ISM: molecules - ISM: structure - 
radio lines: ISM. 


1 INTRODUCTION 


The processes surrounding the life and death of massive 
stars play an important part in the evolution of galax¬ 
ies at all epochs (see e.g. Hennebelle fc CommerconI I2OI2I : 


ISafranek-Shrader et allI2OI4I '). while at all times turbulence 

in the interstellar medium (ISM) pl ays a predomina nt role in 
regulating massive star formation (lFederrathll2013l b Hence, 
understanding how turbulence shapes star formation, and 
how star formation in turn contributes to driving interstel¬ 
lar turbulence, is an important step in understanding diverse 
phenomena such as the evolution of the molecular ISM in 
galaxies, the formation of massive stars, and, eventually, the 
role that turbulence may play in the formation of planetary 
systems. 

Three particularly important topics for understanding 
the interaction between turbulence and massive star for¬ 
mation are: i) The sources of the energy required to drive 
interstellar turbulence; ii) their relative importance at dif¬ 
ferent scales (with large-scale Galactic flows, supernova ex¬ 
plosions, outflows from young, massive stars and expand- 
ing Hii regions all likely to contribute at various s cales e.g. 
ImA Low fc Klessenll20M : lMcKee k, Ostrikerj|2007ll . and iii) 
The effect that turbulence and energy injection may have on 


enhan cing or disrupting star formation at large, spiral-arm 
scales jLuna et al.ll2006l'l and smaller giant-molecu lar-cloud- 
size scales (see e.g. Harper-Clark fc Murravll201ll b 

The study of infall and outflow in massive star form¬ 
ing regions is well connected to the subject area of turbu¬ 
lence. Outflow from massive stars ma y contribute to the 
driving of turbu lence in the ISM (e.g. iRivilla et all l2013l : 
iFederrath et al.ll2014l ~l. while turbulent fragmentation of gas 
that is infalling on to a protostellar cluster may change the 
number and mass distribution of the stars form ing in the 
cluster jPeters et al.ll201ol : iGirichidis et al.ll2012|j . 

Bipolar molecular outflows are foun d ubiquitously 
across all form i ng stellar size scal es (e.g. ISu et all l2004l : 
IwiT et al I I2OO5I: IZ harig et al.l 120131') down even to brown 
dwarfs (IWhelan et al.l 120071 1. The mass of molecular ma¬ 
terial observed in the outflows from high mass star form¬ 
ing regions (HMSFRs) is sufficiently large that it is likely 
entrained from the surrounding interstellar medium (ISM) 
in additio n to the component asso ciated with the form¬ 
ing star (iKlaassen fc Wils^ I2OO8I I. The driving mecha¬ 
nism for outflows from HMSFRs remains unclear but may 
also be due to a similar process to that of low mass star 
formation. In the accretion model of star formation (e.g. 
IShu. Adams &: Lizand[l987l l. gravitational infall of the sur- 
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rounding material onto a disc snrrounding a forming star 
leads to mass ejection and the dissipation of excess angu¬ 
lar momentum in the form of a bipolar jet; the causal¬ 
ity of infall and outflow means whe re one is detected, the 
other is likely to be present. The IShu. Adams fc Lizand 
(Il987ti model considers isolated star formation; modelling 
suggests that outflows associated with high-mass stars can 
be hi ghly collimated. In an interesting and recent develop- 


Mac Low et al.l (l2014h show that when multiple stars 


ment. 

form in a common accretion flow, such as accretion onto a 
massive protostellar cluster, many star and protostars within 
the cluster have common outflow axes. Hence, even with 
massive star formation, there is likely to be a strong associ¬ 
ation between infall and outfl ow. Observational evidenc e for 
this scenario can be found in iKlaassen fc WilsorJ ll2008lj . 

To provide benchmarking observational constraints for 
the role of turbulence in giant molecular clouds, a multi- 
molecular line mapping of the G333 Giant Molecular Cloud 
(GMC) with t he Mopra and Nanten2 telescopes ha s been un¬ 
dertaken (see Bains et al.ll200(j : I Wong et al.ll200^ : ILo et all 


l2007l . [200a . l201l|j . In collecting data for this project, we have 

serendipitously detected signatures of outflow and infall in 
three of the brighter molecular features in the GMC, which 
we have designated G333.6-0.2, G333.1-0.4 and G332.8- 
0.5. Within a Mopra beam size of the outflow sources are 
radio-detected Hll reg ions, 1.2-mm dust emission clumps 
l|Mookeriea et al1l2004r i. CH3OH and HjO masers and IRAS 
sources, features which are all consistent with the presence 
of high mass star formation. 


1.1 G333 outflow/infall sources 


In this paper we examine evidence of infall and outflow to¬ 
wards three massive star forming sources in the G333 mas¬ 
sive molecular cloud complex (see Figure [T]). The G333.6— 
0.2 molecular outflow is associated with the most massive 
of t he three HMSFRs an d harbours a young OB cluster 
(e.g. Fuiivoshi et al.l l2006lj . with a dust mass of 1.6 x 10^ 
M^ llMookeriea et ahl 12004 1. The IRAS source closest to 
G333.6-0.2, IRAS 16183-4958, is one of the most lumi¬ 
nous far-infrared (FIR) sources in the sky. From observa - 
tions with the VLT MIR VISIR camera, I Grave et al.1 (l2014| j 
found indications that this region consists of two main lu¬ 
minous sources (04V and 05V) which account for at least 
half of the luminosity from this region. As well as being 
proximate to the Hll region G333.6-0.22, this molecular out¬ 
flow sou rce is also associated with MSX and 1.2-mm dust 
sources ( |Mookerjea et al] 200j) and masers of HjO and OH 
dBatchelor et al.l 1980l : Caswellll998l j. 


The HMSFR associated with G333.1—0.4 is interme¬ 
diate in mass to the other two sources discussed in this pa¬ 
per. High resolution deep near-IR imaging and spectroscopy 
reveals an emb edded OB star cluste r in very early evolu¬ 
tionary stages (iFigueredo et al.ll2005l j. Spectroscopy shows 
the two brightest stars in the cluster have spectral type of 
06 and 08 stars, and numerous young stellar objects with 
excess near-IR emission due to circumstellar disks or en¬ 
velopes. The cluster has integrated mass of 1.0 x 10® M©. 
This source also has the most prominent SiO emission in the 
whole G333 comple x, aside from th e cold core G333.125- 
0.562 as discussed in lLo et"^ I (l2007l) . 

Finally, we also observed G332.8—0.5 which is the 


sma llest of the three HMSF Rs, with a dust mass of 5.5 x 10® 
Mq llMookeriea et al.ll2003 ~l. It has FIR colour characteris - 
tics of an ultra compact Hll region dBronfman et al.lll996l l. 
Table [ 1 ] gives a summary of previous observational identifi¬ 
cations of the three sources in infrared continuum and maser 
transitions. 

The molecules presented in this paper trace different 
conditions in the interstellar medium related to infall, out¬ 
flow and the dense cores in the sources. Being such a ubiq¬ 
uitous species, the ®^GO emission traces the full spatial and 
kinematical extent of the outflows. The ®®GO traces them 
to a lesser extent, but where present it can be used with the 
®^GO to calculate the molecular column density (under the 
assumption that the ®®CO is optically thin). Emission from 
C®®0 is mainly confined to the cloud cores and provides a 
measure of the core mass, if assumed to be optically thin. 
SiO emission is known to be enhanced in outflows due to the 
presenc e of shocks, although it is n ot found in every outflow 
source dXlaassen fc Wilson 2007 ). HGO~*~ can trace both 
infall and outflows (e.g. Rawli ngs et al.ll2004l : iMvers et al.1 
ll996l : lKlaassen fc Wilsonl2007^ as it easily becomes optically 
thick and traces a high critical density. Emission from the 
SiO and HCO^ lines therefore provides a tracer of outflow 
and infall phenomena. GS is a high density (~ 10®cm“®) 
tracer which is found towards star-forming condensations 
rather than outflow wings and so traces the systemic ve¬ 
locity of the clouds and the degree of turbulence. NjE"*" is 
prominent in cold, dense cloud cores rather than in outflow 
wings, typically at Tex 5, ‘20 K, due to its main destr oyer 
CO being depleted (iBergin. Langer fc Goldsmithlll9^ 1. 

Eor the three G333 sources, we present Mopra data 
of molecular line emission which show evidence for outflow 
and infall, specifically that of the low excitation rotational 
transitions of three CO isotopologues, CS, two isoptopo- 
logues of HCO"*", SiO and NjH^. The observations are de¬ 
scribed in Section [21 In Section O we present archival Spitzer 
GLIMPSE mid-IR imagery of the three outflow sources over¬ 
laid with contours of Mopra CS (J = 2 — 1) data, followed 
by the Mopra spectral line profiles and discuss the evidence 
for the presence of infall from these lines. In Section (4) we 
use the CO isotopologue data to calculate the column den¬ 
sity in the outflows and use this to derive their mass and 
energetics. In Section [5] we summarise the observationally 
derived properties we have dete rmined for these so urces and 
introduce a companion paper (IWiles et al.ll2015h in which 
models of the regions are presented. 


2 OBSERVATIONS AND DATA REDUCTION 

The molecular line data presented here are comprised of 
data from our Mopra G333 multi-molecular lines m apping 
llBains et al.ll2006l : IWong et al.ll2008l : ILo et al.ll2009ll as well 
as ®^CO maps of the individual sources. For a more detailed 
description of the observing procedure and data reduction 
steps, we refer the reader to the references mentioned im¬ 
mediately above. 

The Mopr4i] radio telescope is a 22-metre-single-dish 


^ The Mopra radio telescope is part of the Australia Telescope 
National Facility which is funded by the Commonwealth of Aus¬ 
tralia tor operation as a National Facility managed by CSIRO. 
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Figure 1. Colour-composite images of the Spitzer IRAC photometric bands for the vicinity of the outflow sources where red is S.O-pm, 
green is 4.5-/i.m and blue is 3.6-/i.m. Overlain are white contours of total intensity Mopra CS (J = 2 — 1) emission, which picks out the 
cores of the outflow sources. Contour levels start at 10 K kms~^, with increments of 5 K kms“^ up to 35 K kms“^. (a) Overall image of 
the G333 region with the positions of the three sources in this paper indicated. Zoom in image of the three regions: (b) G333.6-0.2, (c) 
G333.1-0.4 and (d) G332.8—0.5. The crosses mark the positions for the spectra shown in Figure[2] as listed in Table[5] which are also the 
peaks of CO emission. The angular separations between the S-pm and CO peaks are 20, 15 and 7 arcseconds for G333.6-0.2, G333.1—0.4 
and G332.8—0.5 respectively, note that the beam size of the Mopra data is 36 arcseconds, thus all within one bean size. (Colour image 
in the electronic version.) 


Table 1. Names associated with continuum and maser sources found within 80 ar cseconds from the source s fro m the SIMBAD data base. 
The Mopra beam is 36 a rcseconds at 3-mm wav elen gths. The ref e rence s are : fa) iMookeriea et al ] l|2004| ); (b) Becklin et al.l <ll973h : (c) 
iMurphv et al.l l|2QlCtl : fdi iFieueredo et al.l l|2005l '): fel iBreen et al.l l|20Q7lh ffi lCaswelll i|l998l'l : Icaswell et al. ~ i^995l 'l. 


Source 

RA, Dec 
(J2000) 


IRAS 

MSX 

1.2-mm®’ 

Hii 


OH* 

CHgOHS 

G333.6-0.2 

16:22:09.0, 

-50:06:21 

IRAS 

16183-4958 

G333.6046-00.2124 
G333.6044-00.2165 

MMS5 

G333.60-00.21^ 

G333.6080.215 

G333.608-00.215 

- 

G333.1-0.4 

16:21:03.3, 

-50:35:12 

IRAS 

16172-5028 

G333.1104-00.4223 
G333.1306-00.4257 

MMS39 

MMS40 

G333.1-00.4'^ 

G333.1210.434 

G333.135-00.431 

G333.121-00.434 

G332.8-0.5 

16:20:08.9, 

-50:53:14 

IRAS 

16164-5046 

G332.8269-00.5489 

MMS68 

G332.8-00.6‘^ 

G332.8260.549 

G332.824-00.548 

- 
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Table 2. Molecular lines observed. Columns are as follows: (1) molecule, (2) transition; (3) rest frequency (for NjH'*" only the frequency 
of the main hyperfine line is given); (4) months for observations; (5) velocity resolution; (6) typical Icr rms off-line noise level per velocity 
channel in terms of measured T^; (7) reference. 


Molecule 

(1) 

Transition 

(2) 

Rest frequency 
GHz 
(3) 

Months Observed 

(4) 

Av 

kms“^ 

(5) 

ICT 

K 

(6) 

Reference 

(7) 


J 

= 1-0 

115.27 

2006 Aug 

0.09 

0.2 

this work 

13CO 

J 

= 1-0 

110.20 

2004 Jun - Oct 

0.17 

0.1 

Bains et al. ('2006^ 

ci®0 

J 

= 1-0 

109.78 

2005 Jul- Sep 

0.17 

0.1 

Wong et al. 120081 

CS 

J 

= 2-1 

97.981 

2006 Sep - Oct 

0.10 

0.1 

Lo et al. ('2009') 

HCO+ 

J 

= 1-0 

89.190 

2006 Jul- Sep 

0.11 

0.1 

Lo et al. ('2009') 

h13co+ 

J 

= 1-0 

86.754 

2006 Jul - Sep 

0.12 

0.2 

Lo et al. 12009) 

N2H+ 

J 

= 1-0 

93.173 

2006 Sep - Oct 

0.11 

0.1 

Lo et al. 12009) 

SiO 

J 

= 2-1 

86.847 

2006 Jul- Sep 

0.12 

0.1 

Lo et al. 12009) 


telescope located near Coonabarrabran, NSW, Australia, 
with a beam size of 36 arcseconds at 3-mm wavelengths. The 
main beam brightness temperature Tmb and antenna tem¬ 
perature TX are related by the antenna efficiency r]i, at fre¬ 
quency V such that Tmb = and this was us ed to derive 

Tmb . The Mopra beam has been characterised bv iLadd et al.l 
(l2005h and the beam efficiencies used are as listed there. 
The observing bandwidth was conhgured so that the cen¬ 
tral channel corresponded to —50 kms“^, the approximate 
velocity at which the emission from the GMC complex is 
centred. The refer ence (OFF) positi on is at q:j2ooo = 16 : 27, 
5 j2000 = —51 : 30 jBains et al.l200^ . Throughout this work, 
velocities are given in the radio convention and in terms of 
uiar, that is, with respect to the kinematic local standard 
of rest (Isr). The data were reduced using the livedata 
and GRIDZILLA packages available from the CSIRO/CAS^, 
weighted by the relevant Tsys measurements, and have been 
continuum subtracted. We summarise the observational de¬ 
tails in Table (2] The G333 cloud was observed in a number 
of other sp ecies in additio n to those listed in section 1.1, as 
detailed in lLo etld] (l2009l ): however, these are not discussed 
in this paper. 

3 RESULTS 

3.1 Mopra Molecular Line Data 

3.1.1 Velocity Profiles 

In Figure [2l we show the molecular line velocity profiles 
taken at the spatial location of the peak GO emission for 
each source (marked with red crosses in Figure [T] and listed 
in Table [5]). The displayed velocity range of the profiles was 
determined from the data, the line which traces the 

maximum extent of the outflows; it is this velocity range 
that was used to perform the analysis of the outflows as 
described in Section [T] 

In each source, multiple velocity features are visible in 
some line profiles, particularly those of and HCO”'", 

whilst GS, G^®0 and show a single velocity fea¬ 

ture. GS is a high-density tracer and thus associated with 
the denser cloud core region, while the H^^GO’*' isotopo- 
logue is assumed to be optically thin and hence also traces 
the core rather than the outflow. The C^®0 emission also 

^ See URL http://www.atnf.csiro.au/computing/software/livedata 


appears to mainly trace the cloud core component. is 

found in cold dense cores rather than outflows. In the case 
of G333.6-0.2 and G332.8-0.5, we consider these to be asso¬ 
ciated with a single source rather than several overlapping 
velocity features along the line-of-sight. This is evidenced by 
the CS, G^®0 and H^^GO’*’ line profiles, which comprise a 
single component at the same central velocity. 

In both G333.6-0.2 and G332.8-0.5, the and 

HCO^ lines all show the broad wings characteristic of out¬ 
flows, particularly on the blue-shifted side in the case of 
G333.6--0.2 (detailed analysis see Section 14.211 . In G333.1- 
0.4, a blue-shifted shoulder (at ~ —60 kms~^) is present in 
all three CO isotopologues and also HCO’*'; the fact that 
this feature is present in both optically thick (CO) and thin 
(0*^*0) lines, unlike the other two sources, may suggest that 
it is due to a confusing cloud along the same line-of-sight 
rather than part of the wing emission in G333.1--0.4. The 
detection of emission from HCO"^ and SiO i n all the sources 
is consistent with the presence of outflows dRawlings et al.l 
|2 oo 4 [iu aassen fc WilsonllioO^) . 

An infall signature in a spectral line presents itself in 
the form of a red-blue asym metry, usually with a diminished 
red-s hifted component (e.g. lWalker et aPl 19941 : [Mvers et al.l 
Il996l) . Such an asymmetry is apparent to varying degrees 
in the line profiles of all three sources. G333.6-0.2 shows 
the most extreme and broadest red-shifted infall feature, 
extending from ~ —33 kms“^ bluewards towards the centre 
of the line in and HCO”'" (Figure [2J and also visi¬ 

ble to a lesser degree in ^®CO. Continuum absorption has 
further ‘distorted’ the line shape of HCO"^ such that the 
infall feature is negative (the profiles have been continuum 
subtracted). This is consistent with the presence of a num¬ 
ber of radio- and millimetre- wavelength detected c ontinuum 
sources within G333.6-0.2 (iFuiivoshi et al.l l2006l . see also 
Figure [TJ. Indeed, the peak in the continuum absorption 
of the molecular line emission occurs along the line-of-sight 
to the position of peak radio flux density to the Hll region 
G333.6-0.22 (observed with the Australia Telescope Com¬ 
pact Array), which also coincides with the 1.2-mm dust peak 
(Figure O. 

G333.1-0.4 also clearly shows an abrupt fall-off on the 
red-shifted side of its line profiles (Figure [2]). This is a less 
extreme infall feature than that seen in G333.6-0.2 and is 
again consistent with the line being absorbed by the contin¬ 
uum source which we have detected at radio and millimetre 
wavelengths (see also Figure [T]). G332.8-0.5 shows a more 
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Figure 2. Velocity profiles of the molecular lines for each of the outflow sources. The spectra were taken from the positions indicated by 
the red crosses in Figure^ (as listed in Table [ 5 J 1 . Left: G333.6-0.2; centre: G333.1-0.4; right: G332.8-0.5. From top to bottom the lines 
are (see Table[2]l ^^GO, ^^CO, G^®0, GS, HCO”^, N 2 H'*' and SiO. The y-axis is the intensity, T^, in units of K and the x-axis 

is the v\sr velocity in kms“^. The dashed lines indicate = 0 (horizontal) and centroid velocity (vertical). 
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Figure 3. Total intensity of HCO+ greyscale (Mopra Tele¬ 
scope) with dashed contours, and black contours overlaid from 
the G333.6—0.22 Hll region at 4.8-GHz (ATCA). The lower left 
circles indicate the beam sizes of the data, 36 and 6 arcseconds 
for Mopra and ATCA respectively. 

‘classic’ infall profile, with a clear split of ~ 5 kms“^ be¬ 
tween asymmetric red- and blue-shifted peaks (—54 and —60 
kms~^ respectively) which is particularly evident in the line 
profiles of and HCO+. 

3.1.2 Outflow Maps 

In FigureUl we present images of the data to show the 

overall velocity structure. The zeroth moment (i.e. total in¬ 
tensity) contours of the red and blue wing emission (summed 
over the velocity ranges defined in Table 0 for each out¬ 
flow source are shown overlaid with symbols indicating the 
positions of other likely related sources of emission in the 
region, more details of which are given in Table [T] Two of 
the three sources (G333.2-0.6 and G333.1-0.4) have a clear 
offset between the red- and blue-shifted total intensity emis¬ 
sion; the position angle (PA) of this offset is different for each 
source, i.e. there is no evidence for a general alignment of 
such kinematic characteristics across the large-scale G333 
GMC. However, care must be taken in the interpretation 
of the red-shifted total intensity emission from G333.1-0.4 
and G333.6-0.2 because of the red-shifted absorption dip 
in each of these sources. For source G332.8-0.5, there is no 
measurable blue-shifted emission thus it is not possible to 
determine the outflow axis. 

We also show the PV (position velocity)-arrays in the 
right column of Figure [4] , made by taking a slice through 
the ^^GO data cube along a plane defined by the positions of 
the peak emission in the most extreme red- and blue-shifted 
channels; the slice positions are indicated by the thick solid 
lines overlaid on the moment images. The PV-arrays for 
G333.6-0.2 and G333.1-0.4 show the outflows blending into 
spatially and/or kinematically adjacent emission, and the 
abrupt falls in the red-shifted channel emission. The G332.8- 
0.5 PV-array is not shown because the outflow axis is not 
determined. 

3.2 Spitzer GLIMPSE Imagery 

In Figure[Tl we show the Spitzer Space Telescope GLIMPSE 
(Galactic Legacy Infrared Mid-Plane Survey Extraordinaire; 


Table 3. Gaussian fits to the C^®0, and GS line profiles 

of the cloud core emission in each source (1). Both the centroid 
(2), (3), (4) and FWHM (5), (6), (7) velocities are given for each 
line. 


Source 


(km s 


FWHM (km s- 



ci®0 

13G0 

CS 

ci®o 

13CO 

CS 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

G333.6-0.2 

-47.7 

-48.6 

-47.7 

6.1 

7.5 

7.0 

G333.1-0.4 

-52.7 

-52.6 

-52.2 

6.1 

7.4 

6.8 

G332.8-0.5 

-58.5 

-58.3 

-57.7 

5.9 

7.7 

7.0 


iBeniamin et al.l l2003f) 3-colour composite infrared images 
of the whole cloud complex and the three bipolar outflow 
sources. The images are overlain with contours of the CS 
total intensity, which we display here as it has a high as¬ 
sociated critical density (~ 10®cm“®) and so picks out the 
structure of the high-density molecular cores well. 

Infrared emission associated with all three outflow 
sources is apparent in the GLIMPSE images. G333.6-0.2 and 
G333.1-0.4 display complex structures comprising a number 
of knots and filaments in the IR while G332.8-0.5 has a sim¬ 
pler, bipolar nebulosity. The predominance of 8-/mi emission 
(red) in all three sources is consistent with the presence of 
polycyclic aromatic hydrocarbons (PAHs) typically found in 
photodissociation regions (PDRs) around HMSFRs, so the 
8-^m peak is taken as being a sign of an Hll region. The po¬ 
sitions of the Hll regions visible in the Spitzer data are con- 
sistent with the posit ions of the radio continuum emission 
dUrauhart et aLll2007l) . Emission a ssociated with the m as- 
sive OB cluster in G333.6-0.2 (e.g. iFuiivoshi et al.|[2006l ) is 
clearly visible as the saturated region in the Spitzer InfraRed 
Array Camera (IRAC) image and is clearly offset from the 
molecular emission peak, which is itself associated with a 
dark region in the IR emission. Similarly, the molecular peak 
in G333.1-0.4 is offset from the emission associated with the 
Hll region and coincident with a trough in the IR emission. 
Such offsets between molecular emission and Hll regions are 
consistent with more evolved HMSFRs whose feedback ef¬ 
fects have cleared their natal environs. Conversely, the Hll 
region in G332.8-0.5 is located in the centre of the peak 
molecular contour, suggesting it may be a younger source 
and/or less powerful. 


4 ANALYSIS 

4.1 Molecular Lines: optical depths 

As the three HMSFRs are located within a massive GMC 
complex, the analysis was hampered by the multiplicity of 
sources present and their associated confusion, both spa¬ 
tially and kinematically. The outflows are thus less well- 
defined than in low mass star forming regions (LMSFR) and 
the clean, detailed analysis that is possible there is not pos¬ 
sible here. 

We defined the systemic velocity of the outflow sources 
and their core/outflow velocity boundaries by fitting 1- 
dimensional Gaussians to the C^®0, ^®CO and CS line pro¬ 
files taken at the positions along the line-of-sight to the out¬ 
flow centres. These lines appear to trace the core features 
only and have little emission in the wings. The systemic and 
FWHM velocities so determined are listed in Table [3] It can 
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RA (J2000) 


Figure 4. (left:) Total intensity red- (dashed contours) and blue-shifted (solid) wing emission, summed over the velocity ranges 

given in TableFor G333.6-0.2, contours are plotted in 5 K kms~^ steps starting at 5 K kms~^ for both blue and red shifted outflow 
emissions. For G333.1-0.4, contours are plotted in 2 K kms~^ steps starting at 14 K kms“^ for the blue-shifted emission, and in 1 
K kms“^ steps starting at 8 K kms“^ for the red-shifted emission. For G332.8—0.5, there is no detectable blue-shifted emission, the 
red-shifted emission contours are in steps of 2 K kms“^ starting at 18 K kms“^. The thick solid line on each plot indicates the position 
of the slice taken to produce the PV-arrays shown, (right:) Position-velocity arrays with contours plotted at 10 per cent steps of the peak 
emission, (top:) G333.6-0.2. PV-array centred on RA 16:22:9.2, Dec -50:06:02, PA 34°. (middle:) G333.1-0.4. PV-array centred on RA 
16:21:0.6, Dec -50:35:16, PA 351°. (bottom:) G332.8-0.5. PV-array is not shown as the PA is unknown. All angles are measured E from 
N. The symbols indicate the other sources of emission in the region (see Tableas follows: MSX-hox, H 2 O maser-filled triangle, OH 
maser-A, CH 3 OH maser-*, Hii region-©, 1.2-mm source-* and the IRAS sources are indicated by their error ellipses. 


© RAS. MNRAS OOO.flM^ 



















8 N. Lo et al. 


Table 4. Line-derived excitation temperatures, isotopologue intensity ratios and optical depths. The excitation temperature Tex (2) is 
determined from the brightest value of found in each source (1), divided by at the visr velocity listed (3). The maximum 

value for the isotopologue ratios ^^CO/^^CO (4) and ^^CO/C^^O (6) are also listed, together with the corresponding uisr velocities (5), 
(7) where this occurs. The maximum value for the optical depth (8) is determined from the minimum value of the ^^C0/C^®0 ratio, 
and occurs at the velocity (9) indicated. 


Source 

(1) 

Tex 

K 

(2) 

«^lsr 

km 

(3) 

l^GO/l^GO,nax 

(4) 

'i'lar 

km 

(5) 

l^CO/C"»Omax 

(6) 

‘^Isr 

km 

(7) 

'^max 

(8) 

km 

(9) 

G333. 6 - 0.2 

30 

-50.4 

6.4 

-55.7 

25.2 

-52.2 

73 

-44.7 

G333.1-0.4 

33 

-50.9 

4.4 

-62.9 

24.8 

-57.4 

20 

-53.6 

G332.8-0.5 

17 

-59.7 

4.3 

-51.6 

13.8 

-53.3 

37 

-58.3 


Table 5. Measured outflow parameters. The columns are as follows: (1) source; (2) Sz (3) RA and Dec (J2000) giving the position of 
spectra shown in Figure [2] (4) Sz (5) Sz giving velocity limits for blue-shifted outflow; (6) Sz (7) velocity limits for 

red-shifted outflows; (8) Sz (9) 'WBinean ^ ^R-mean mean projected velocity of blue- and red-shifted outflows; (10) position angle for 
the outflow in the plane of the sky, measured E from N. Outflow velocities are determined after subtracting the scaled Gaussian fits to 
the line from the extinction-corrected line profile (see text) (note: no blue-shifted emission is detectable in the line profile 

for G332.8—0.5). The line core emission is taken as between v-q^ and . Parameters calculated for the core and outflow given in other 
Tables use these velocity limits, together with the line centre velocities given by the fit to the profile (see Table 0. 


Source 

RA 

Dec 


UB2 

■"Ri 

«R2 

^Bmean 

^Rmean 

Position 


J2000 

km 

km 

km s ^ 

km 

kms“^ 

km 

Angle (°E of N) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

G333. 6 - 0.2 

16:22:09.0 

-50:06:21 

-69 

-57 

-37 

-26 

15.3 

16.2 

34 

G333.1-0.4 

16:21:03.3 

-50:35:12 

-75 

-63 

-43 

-37 

16.3 

12.7 

351 

G332.8-0.5 

16:20:08.9 

-50:53:14 

- 

-68 

-49 

-37 

- 

15.5 

- 


Table 6. Line fluxes for outflows and core, and their Itr errors, 
Tmb in Kkms“^ (i.e. corrected for beam efficiency), for the 
three CO isotopologues measured in each source (1). For the 
line the fluxes for the blue (2) and red-shifted (4) outflows are 
determined over the velocity ranges defined in Table [5] The line 
core velocity range used for the fluxes of (3), (5) and 

C^®0 (6) lines is between ubj a-nd ub^ in Table [5] For each flux 
the corresponding error is listed. The errors include the random 
errors and an estimate in the uncertainty for the continuum level. 


Source 

(1) 

Blue 

(2) 

Core 

(3) 

Red 

(4) 

Core 

(5) 

C^'^O 

Core 

(6) 

G333.6-0.2 

41 ± 3 

218 ± 4 

42 ± 2 

66 ± 3 

15 ± 2 

G333.1-0.4 

42 ± 3 

369 ± 5 

15 ± 2 

102 ± 3 

22 ± 2 

G332.8-0.5 

- ± 2 

174 ± 4 

39 ± 3 

68 ± 3 

14 ± 2 


Table 7. Fluxes for non-CO lines in Table [2] together with their 
Icr errors, Tmb in Kkms“^ (i.e. corrected for beam effi¬ 
ciency), measured in the core component of each source. The er¬ 
rors include the random errors and an estimate in the uncertainty 
for the continuum level. The velocity range used is between UB 2 
and urj in Table [S] 


Source 

CS 

HCO+ 

H^^CO+ 

SiO 

N 2 H+ 

G333.6-0.2 

32±2 

19±2 

6.3±3.1 

1.7±1.3 

8.3±1.3 

G333.1-0.4 

53±2 

35±1 

8.4±2.5 

4.2±1.3 

19±3.7 

G332.8-0.5 

24±2 

15±1 

4.2±1.5 

1.3±1.1 

14±1.4 


be seen that for each source the Usys measurements agree to 
within 1 kms“^. In addition, the and CS lines are seen 
to be moderately optically thick, evident by the ~ 1 kms“^ 
larger values for their FWHM. 

The CS emission in the core of G333.1-0.4 (Figure [2|) 
appears to be comprised of two blended velocity components 
which could not be fit with a single Gaussian. The C^®0 
emission has a smoother core profile and provided a better fit 


in this source. For consistency in our analysis, we therefore 
use the parameters determined from the fits to the C^®0 
profiles in the calculations below. We measured the outflow 
masses and energetics using the observed brightness ratios 
between the *^^00, ^®CO and 0*^*0 line emission. Where 
these lines are well detected (S/N > 5) we may use the ratio 
with the main line to determine the optical depth of the 
emission. We summarise here our use of a radiative transfer 
analysis of the observed emission. 

The radiative transfer equation has the general solu¬ 
tion for the intensity of a velocity cha nnel in a spectral line , 
known as the Detection Equation (e.g. lStahler fc Pallall2005l . 
Appendix C): 

T = TX/ 77 . = /[J.(rex) - MTBom - e-^n (1) 

where TX is the measure d intensity and rji , the beam effi¬ 
ciency (taken as 0.55 from iLadd et al.l (l2005l )). / is the beam 
filling factor for the emission, Jv{T) = [hu — 1] 
with Tex being the excitation temperature and Tbg being 
the temperature of the cosmic background radiation (i.e. 
2.726 K). Tv is the optical depth of the emission. 

Generally, in using Equation [T] we assume for the fine 
isotopologue pair, say *^^00 and ^®CO, that the former 
is optically thick, while the latter is optically thin every 
where. However unlikely to be optically thin in the 

cores of these sources, and thus the derived core masses 
are a lower limit. We also assume they have the same ex¬ 
citation temperature and beam filling factor. Eurthermore, 
ri2co/Ti3co = V(i2co)/Ar(i3cO) = X[^'^C/^^C\, the 
abundance ratio of to ^®CO. We have adopted abun¬ 

dance ratios of X\^^C/^^G\ = 20, as preli minary generic 
radiat ive transfer models undertaken as per ICarolan et al.l 
(l2008l) did not support a higher ratio, if the observed fine 
intensities were to be reproduced. The follow-up modelling 
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tailore d to each of the so urces is discussed in the companion 
paper dWiles et alj 1201^1 to this one mentioned in the In¬ 
troduction (Section [T]). For the oxygen isotopes, we adopt 
a ratio of = 100. These values are smaller 

than the broad Galactic scale values dMilam et al.l[2005ll . 
but taken i nto account of fractionatio n observed in cold 
dens e core dMla denovic fc Roue^ 20W ) an d are c onsistent 
with ICarolan et al. ( 2008l l aiid^ Wiles et al.l d2015h . Assum¬ 
ing these values, we see from the Detection Equation that 
— TA(^^CO)/rX(’^^CO), with T12CO then determined 
by the preceding formula (and similarly for tqisq). 

We hence apply this analysis to determine the opti¬ 
cal depth of the three CO isotopologues lines, as a func¬ 
tion of velocity. Furthermor e, from standard molecular ra¬ 
diative transfer theory (e.g. IColdsmith fc Langerl [l9^ 'l we 
may show that the column density of the upper level of each 
transition is given by 


Nu 


T 

Ahc^ 


An 


( 2 ) 


where k, h and c are the well-known physical constants, v 
is the frequency of the transition, A is the radiative decay 
rate, and An is the channel velocity spacing. The optically 
thin case (r << 1) simply has the optical depth correction 
factor r/(l — e“’^) set to unity. 

Once Nu has been determined the total molecu¬ 
lar column density can be found by applying Ntot = 
(Au/ 5 u)Q(T'ex)e^“/^^®* where = 2J + 1 — 3 is the level 
degeneracy and the partition function is QiTeyl) = 2rex/Tb 
(with To = hvik) at an assumed excitation temperature, 
Tex. We take the values of Tex given in Table IH which are 
the peak values of TX/r]v measured for the ^^CO line in 
each source. The observed ^^CO lines are extremely opti¬ 
cally thick, and it is possible that the actual Tex is higher 
than the temperature derived from the ^^CO flux, thus un¬ 
derestimating the total molecular column density. For a Tex 
of 60 K, the total column density would be almost double 
that of an excitation temperature of 30 K. 

In TableUwe also list the maximum value of the optical 
depth found in each source, which corresponds to the mini¬ 
mum value of the ^^CO/C^®0 ratio, and the corresponding 
velocity. We also list the maximum values of ^^CO/^^CO and 
^^CO/C^®0 found for each source and their corresponding 
velocity. These are seen to be very much less than the as¬ 
sumed abundance ratio, supporting the assumptions made 
above regarding the determination of the optical deptI0. 


4.2 Molecular Lines: outflows, infall and 
turbulence 

In order to determine the velocity extent of the outflow, we 
use the optical depth of (as calculated in Equation 

[T]) to plot a core Gaussian line shape from G^®0, and then 
from it to identify the outflow channels. The remaining pro¬ 
file (see Figure [5]) shows the presence of extended red- and 
blue-wings in each source (except G332.8-0.5 for which no 
blue wing is evident). We use these profiles to determine the 


® Note also that any parameters derived from the line ratios are 
only calculated when the S/N in both corresponding lines is > 5. 



Velocity (km/s) 


Figure 5. Optical depth-corrected, subtracted profiles, showing 
the outflow wings for each source. From top to bottom are shown 
G333.6-0.2, G333.1-0.4 and G332.8-0.5, respectively. The dotted 
line shows the line profile and the dashed line the Gaussian 

fit to the profile, respectively, scaled to the peak intensity 

of the optical-depth corrected line. The solid line shows 

the optical-depth corrected profile with this scaled fit sub¬ 

tracted off (note that in the line core, where this subtraction is 
imperfect due to the scaling, it has been blanked out). This line 
represents the outflow profile. The y-axis scale is in K and the 
x-axis is the velocity in kms“^. 


velocity extent of the outflows, as listed in Table (5] as well 
as the line core, for each source. Table [6] lists the CO line 
fluxes that have then been determined over these velocity 
ranges, corrected for the beam efficiency. Table [3 provides 
the fluxes for the other fines observed in the line core ve¬ 
locity range. In Table [8] the corresponding column densities 
are listed, calculated according to method described above 
using the CO isotopologues. The Table includes the upper 
state column density, Nu, for the *^^00 (i.e. J = 1) and the 
total column density, Atotai- Atotai includes a correction 
for the optical depth listed in the Table (based on the C^®0 
fine) for the integrated flux in the fine core, but assumes 
optically thin emission for the outflow wings (the ^^^CO and 
C^®0 lines are not detected with sufficient S/N in the out¬ 
flow wings to determine their optical depth). A ^^CO/Hj 
abundance of 1.5 x 10“® is also assumed. 

The molecular mass is given by Mh 2 = AtotaiUd^/ruiHs 
where 11 is the solid angle (for which the Mopra beamsize 
of 36” is used), the source distance d = 3600pc, muj is the 
mass of a hydrogen molecule and ^ = 1.2 is the factor taken 
as the mean mass per hydrogen molecule (for fully molecu¬ 
lar gas and helium mass fraction 10 per cent). We estimate 
that the ubiquitous low-level ambient emission may add an 
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error of up to 10 per cent to the measured flux densities and 
hence the masses. A larger possible source of error lies in 
the probable blending of features that is impossible to dis¬ 
entangle without higher excitation transitions and/or higher 
angular resolution. 

We measure a mean length scale I for the outflows 
by calculating the offset of the peak pixels in the extreme 
velocity channels from the centre position of the outflow. 
We then derive a mean outflow time scale foutflow given by 
foutflow = 21/{vb -bUR). Finally, we calculated the quantities 
that characterise the direct mechanical feedback effects of 
the outflows, specifically mass loss rate Mioss, momentum 
p, mechanical force F/iech, mechanical power Lmech, kinetic 
energy Ek and free-fall time tfs: 


Mioss = 


Mb + Mb. 

^outflow 


p = |MbUb| + |MrUr| 


Fmech — 


P 


^outflow 


Ek = -^{MbVb + MrUb) 


(3) 

(4) 

(5) 

( 6 ) 


the MSX and SIMBA datE0. The application of the far-IR 
balloon data is particularly important in determining the 
source luminosity. The fitting applied a 2-component grey- 
body of the form 


Fv = f2hotBi/(2hot) + war m (9) 

where the dust emissivity is given by = 1 — , with the 

optical depth Ti, = to{v/ vo)^, for a dust emissivity index jS 
taken to be equal to 2. corresponds to the wavelength Aq 
where the IR emission becomes optically thin (i.e. = 1). 

Fitting provides estimate for Thot and Twarm, representa¬ 
tive temperatures for the hot and warm components of the 
fit, although only the warm component can be interpreted 
as physical par ameter characteristic of the source (see e.g. 
iHill et al]l2()09l l. is the Planck black-body function. The 
angular sizes, flhot and flwarm, provide an effective source 
size for the IR emission, rhot and rwarm, at the distance to 
the source of 3.6 kpc. For reference, 0.1 pc corresponds to 
an angular size of 6 arcseconds, unresolvable with this data. 
The best fit parameters are listed in Table [TT] 

The fitting also provides a source luminosity (the area 
under the spectral energy distribution), and a dust mass. 
The latter is derived from the optically thin 1.2-mm (250- 
GHz) emission in conjunction with Twarm, the dust temper¬ 
ature determined for the extended, warm component (which 
dominates the total flux) as follows: 


Tmech — 


Ek 

toutflow 


(7) 


tff = ^Ji-Kl'NlGp (8) 

tg requires the average density p, determined from the 
average number density n calculated as described in the next 
section fSection l4.3ll . The results of these calculations for the 
outflow parameters are given in Tabled together with the 
dust-derived mass and luminosities (also taken from Section 

[Oil . 

From Gaussian fits to the C^®0 and CS line profiles, 
the cores have a velocity (ufwhm) of 6 - 7 kms“^ (Table[3|). 
Assuming a gas temperature of 20 K, the Mach number is 
above 12, while for 100 K gas, which is the warm component 
of dust temperature from SED fits (Table Hill , the Mach 
number is 6. Thus, the turbulence of the cores are highly 
supersonic. 


4.3 Dust Continuum 


Spectral energy distributions were fitted to the fluxes for the 
IR emission associated with the three sources (G333.6-0.2, 
G333.1-0.4 and G332.8-0.5). These fluxes were determined 
from a combination of MSX (8 - 21 pm), IRAS (25 - 100 
pm), the Tata Institute of Fundamental Research (TIFR) 
balloon-borne telesc ope measurements at 150 and 210 pm 
llKarnik et aI.l[2001^ and 1 .2-mm emi ssion measured using 
SEST/SIMBA ( Mookeriea et al1l2004h . as listed in Table fTOl 
While the IRAS and balloon measurements used large aper¬ 
tures (3 arcminute for the balloon), inspection of the images 
across these wavebands showed they were dominated by a 
single source, at least within the 20 arcseconds resolution of 


^(dust) = -r (10) 

The total mass opacity coefficient was taken as 
K 250 GHZ = 0.005 g“^cm^ and the gasidust mass ratio as¬ 
sumed to be 100. The core masses derived from the CO 
lines and the dust mass derived from the continuum are 
both given in Table [H] As can be seen, the masses derived 
using the two methods are in good agreement given the as¬ 
sumptions made, differing by less than a factor of 3. 

Luminosities are found to be 5 - 10 x 10® L 0 , dust 
masses a few thousand solar masses, and dust temperatures 
from 70 - 100 K for the three sources. Note that including the 
extended emission around G333.6-0.2 nearly doubles the de¬ 
termined luminosity and dust mass from the region, but the 
other parameters determined are little changed. The dust 
luminosity to mass ratio is found to be ~ 400Lq/Mq in 
all cases. Given the fitted source size it is also possible to 
calculate the average density, column density and optical 
depth for each source, and these figures are also included in 
Table [TT] M(dust), Lfir and tg (derived from n) are also 
included in Table [9| Note that the column density, N in 
Table is applicable to the size of the warm component, 
rwarm. This ueeds to be scaled to the 36 arcseconds Mopra 
beam size to be directly comparable to A/otai in Table (HJ 
which requires multiplying N by ~ 0.02, yielding values 
within a factor of 4 for the column density. 


Note however that for G333.6—0.2 the emission is somewhat 
more extended, so we also repeated the fitting for an aperture 
which includes all the emission within a 3 arcminute aperture. 
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Table 8. Derived column densities and masses. Upper level column density (^^CO J = 1) (3), (4) & (5), total H 2 column density (6), 
(7) h (8) (assuming a [^^CO/H 2 ] abundance of 1.5 x 10“^) and H 2 mass in the blue and red-shifted outflows and in the core component 
(9), (10) Sz, (11), for the three sources (1). For the core component the mean optical depth (2), determined from the ^^C0/C^®0 flux 
ratio, is used to correct for extinction in the determination of A^total mass. The values listed for the outflow components (and 

Supper for tflo core), do not have any extinction correction applied; i.e. optically thin emission is assumed for the outflow parameters. 
The upper limits listed for the blue outflow for G332.8—0.5 assume the Icr error for the flux given in Tabled 

Source Mean r A^upper (10^*^ x cm“^) A^total (lO'^^x cm“^) Mass (Mq) 


(1) 

Line Core 
(2) 

Blue 

(3) 

Core 

(4) 

Red 

(5) 

Blue 

(6) 

Core 

(7) 

Red 

(8) 

Blue 

(9) 

Core 

(10) 

Red 

(11) 

G333.6-0.2 

6.6 

2.7 

15 

2.8 

7.8 

280 

8.0 

34 

1200 

34 

G333.1-0.4 

5.9 

2.8 

25 

1.0 

8.9 

460 

3.1 

38 

1970 

14 

G332.8-0.5 

8.3 

< 0.12 

12 

2.7 

< 0.23 

190 

5.0 

< 1 

810 

22 


Table 9. Derived outflow parameters. The columns are as follows: (1) source; (2) (3) Mb & Mr are the blue and red outflow masses 

respectively; (4) Mcore molecular mass derived from the line core emission; (5) / is the outflow length scale; (6) foutfiow outflow time for 
this scale size, given the outflow speeds listed in Table [S] (7) fff the free-fall time; (8) Mioss mass loss rate derived for each outflow 
(red -h blue components); (9) p total outflow momentum; (10) -Fmech th® mechanical force for the outflows; (11) Tmech their mechanical 
luminosity and (12) F'k their kinetic energy. 


Source 

(1) 

Mb 

Mq 

(2) 

Mr 

Mg 

(3) 

Mcore 

Mg 

(4) 

1 

pc 

(5) 

^outflow 
103 yrs 
(6) 

iff 

10^ yrs 

(7) 

•^loss 

10-3 Mg/yr 
(8) 

V 

MQkms”^ 

(9) 

-^mech 

MQkms”^ yr~^ 
(10) 

T'mech 

L© 

(11) 

Ek 

erg 

(12) 

G333.6-0.2 

34 

34 

1200 

0.17 

10 

6.2 

6.4 

1100 

0.10 

130 

1.7 

G333.1-0.4 

38 

14 

2000 

0.50 

40 

13 

1.5 

800 

0.023 

29 

1.2 

G332.8-0.5 

< 1 

22 

810 

0.67* 

80* 

18 

0.27* 

340 

0.0040* 

5.2* 

0.53 


* since there is no blue outflow detected, the outflow length scale for this source is from the peak red outflow to the core centre, and this number is used for 

subsequent calculation. 


4.4 Luminosities and ionizing photons 

iFuiivoshi et alj ll2006li used hydrogen radio recombination 
lines and the Bracket 7 IR line to derive the electron temper¬ 
ature and so calculate the number of Lyman continuum pho¬ 
tons in G333.6-0.2. They estimate this to be ~ 9.5x 10'^® s“^, 
coming from the inner ~ 50 arcseconds, equivalent to ~ 19 
07V stars from their estimation of the number of Lyman 
conti nuum photons, assu ming ~ 5 x s“^ for one 07V 
star jMajtins_et_aJj 200^1 . 

IConti fc Crowtherj (1200^ 1 estimated an IR luminosity of 
1.1 X 10® LQwith an aperture of 4 arcminute on the MSX and 
IRAS images, compare to 1.1 x 10® LQfrom our SED fittings 
with MSX, TIER balloon and SIMBA data, and 2.0 x 10® 
Lp) if we take a 3 arcmi nute aperture (Table lllll . Both of 
IConti fc Crowtlierl ll2004h and our estimations of luminosity 
are near double the lumino sity from 04V and 05V stars 
found bv jGrave et al.l (l2014ll . consistent with the finding of 
iKumaij (I 2 OI 3 I) that these two stars account for at least half 
of the luminosity of the UCHll region. 


5 SUMMARY 

We have used extensive molecular line data obtained with 
the Mopra radio telescope to search for evidence of out¬ 
flow and infall associated with massive star formation in the 
G333 giant molecular cloud complex. The complexity of such 
sources where widespread massive star formation is under 
way makes such searches difficult. Evidence is generally re¬ 
quired from a variety of molecular species and their isotopo- 
logues to unravel the competing effects of multiple sources, 
optical depth and different evolutionary states within the 
region. We have used data on the 3mm-band emission from 
eight molecular species, including three isotopologues of CO, 
dense gas tracers such as CS and NjH'*", and shock tracers 
such as HCO’*’ and SiO, together with archival continuum 


data from four IR surveys (Spiteer/GLIMPSE, MSX, IRAS, 
TIER balloon), for this purpose. The line data has 0.6 ar¬ 
cminute and 0.1 kms“^ spatial and spectral resolution, and 
we have used it to determine physical parameters for sources 
in G333 to characterise their properties. 

We have identihed three massive star forming sources 
within G333 showing evidence of both infall and outflow: 
G333.6-0.2, G333.1-0.4 and G332.8-0.5. Outflow is evident 
by the broad wings to some of the line profiles, and infall by 
line splitting. These three sources are at different evolution¬ 
ary states, with G332.8-0.5 at the earliest stage of star for¬ 
mation with the IR sources and Hll regions least prominent, 
and vice-versa for G333.6-0.2. G333.1-0.4 lies in-between. 
All sources show broad profiles characteristics of outflows. 
This is particularly prominent in the blue wings for G333.6- 
0.2 and G331.1-0.4, but only evident in the red-wing for 
G332.8-0.5; the blue-wing is absent in this latter source. A 
clear outflow axis can be defined between the offset red and 
blue lobes for the first two sources, though the outflow itself 
remains morphologically poorly defined, in contrast to out¬ 
flows seen in typical low mass star forming regions. HCO”'" 
and SiO emission, prominent in shocks, is also detected in 
all three sources, consistent with outflows being present in 
them. 

Infall signatures are also apparent in the form of a red- 
blue asymmetry in all three sources, the line being self- 
absorbed by a cold, central continuum source. G333.6-0.2 
is the most extreme with a splitting of ~ 15 kms“^, and 
G332.8-0.5 the least at ~ 5 kms“^. 

We have used the ratio of the three GO isotopologues 
to correct for optical depth at each velocity channel, and 
so determine the column density for the line core as well as 
the outflow lobes. From this we are able to determine their 
masses, as well as estimate mass loss rates, outflow mechan¬ 
ical energies and luminosities. Typical outflow masses are 
10 to 40 M 0 in each lobe, compared to core masses of order 
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Table 10. IR source fluxes in Jy measured with MSX^ IRAS^ TIFR balloon-borne telescope llKarnik et al.|[200lf) and SEST/SIMBA 
jMookeriea et al]|20Q4h . The second row for G333.6-0.2 is for fluxes taken within a 3^ aperture. This source is too confused to determine 
reliable IRAS fluxes in its large beam, so these are not included in the fits. Columns (1) and (2) provide our and Karnik’s source ID; (3) 
&: (4) the source position; (5)—(8) the MSX fluxes; (9) - (11) the IRAS fluxes; (12) &; (13) the Balloon fluxes; (14) the SEST/SIMBA 
flux. 


IR Source 

Position 


MSA-band 



JiJ A 5-band 

Balloon 

SIMBA 

Our ID 

Karnik ID 

RA 

Dec 

8.3 

12.1 

14.6 

21.3 

25 

60 

100 

150 

210 

1200 



J2000 











(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

G333.6-0.2 

S23 

16:22:10 

-50:06:00 

461 

3310 

2970 

6810 




27800 

11600 

82 

(3') 




919 

4479 

4223 

12050 




27800 

11600 

128 

G333.1-0.4 

S15 

16:21:01 

-50:35:15 

133 

378 

512 

1550 

3040 

12900 

18400 

20800 

10000 

34 

G332.8-0.5 

Sll 

16:20:11 

-50:53:20 

39 

115 

197 

831 

2030 

14100 

16800 

15400 

7790 

20 


Table 11. Derived source parameters using the fits to equation|9] The second row for G333.6—0.2 is for the fluxes measured through a 3' 
aperture. Columns are as follows: (1) source ID; (2) & (3) dust temperature for the hot and warm components; (4) col d temperature to 
which the temperature falls in the envelope at 1 pc; was estimated from a ID radiative transfer model in DUSTY lllvezic &: Elitzun 

Il997^ (5) &: (6) radii for the hot and warm components; (7) wavelength where the optical depth is unity; (8) source luminosity; (9) dust 
mass (converted into total gas mass); (10) average density; (11) column density (see comment in text for comparison to Table|8}; (12) 
visual extinction; (13) luminosity / mass ratio. 


ID 

(1) 

That 

K 

(2) 

Twarm 

K 

(3) 

Tcold 

K 

(4) 

^hot 
10 “® pc 

(5) 

^ warm 

pc 

(6) 

Ao 

fim 

(7) 

Tfir 

IO^Lq 

(8) 

M(dust) 

M© 

(9) 

n 

10 ® cm-® 
(10) 

N 

lO®'* cm-2 
(11) 

mags 

(12) 

L/M 

LqMq-i 

(13) 

G333.6-0.2 

350 

98 

20 

2.1 

0.081 

752 

1.2 

3300 

29 

20 

9800 

370 

(3’) 

350 

98 

20 

3.3 

0.10 

752 

2.0 

5200 

23 

20 

9800 

390 

G333.1-0.4 

240 

69 

20 

3.3 

0.11 

371 

0.57 

1800 

6.1 

5.6 

2800 

310 

G332.8-0.5 

238 

68 

20 

1.8 

0.11 

285 

0.50 

1100 

3.5 

3.3 

1700 

460 


10® Mq. Outflow size scales are a few tenths of a parsec, 
timescales are several x 10 ^ years and mass loss rates a few 
xl0“^ M 0 /yr. Flow momenta are ~ 1000 MQkms“®and 
their mechanical luminosities a few x 10 Lq . 

The source SEDs were used to calculate their luminosi¬ 
ties, then by fitting to a 2 -component grey-body model also 
the dust mass, dust temperature and source size for the ex¬ 
tended component. Luminosities are ~ 10 ® Lq, dust masses 
a few X 10® Mq (similar to that inferred from the line emis¬ 
sion), dust temperatures ~ 100 K and sizes 0.1 pc. This 
yields number densities of a few x 10 ® cm“® and luminos¬ 
ity/mass ratios, L/M ~ 400 Lq/Mq. The dust luminosity is 
also similar to that inferred from the hydrogen Bry flux for 
G333.6-0.2. 


Parameters for the infall may also be inferred from the 
line splitting in the profiles. However this is a complex pro¬ 
cedure, as the radiative transfer needs to be considered. Es¬ 
timates cannot simply be made from the magnitude of the 
line splitting, but must consider the medium through which 
the radiation passes. This requires a more sophisticated ap¬ 
proach than that presented here, and involves modelling 
of the source geometry as well as its physical characteris- 
tics in order to y ield line profiles. In a companion paper 
dWiles et al.l l2015l ) we apply such a 3D radiative transfer 
analysis making use of the code MOLLIE, which is able to 
consider the competing contributions of the outflow, infall 
and ambient gas, which may also have different densities, 
temperatures and chemical compositions, in order to pro¬ 
vide an estimate of the source parameters, and in particular 
to yield mass infall rates and infall speeds from the data set. 
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